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Abstract. A covariant approach towards a theory of deformations is developed 
to examine both the first and second variation of the Helfrich-Canham 
Hamiltonian — quadratic in extrinsic curvature — which describes fluid vesicles 
at mesoscopic scales. Deformations are decomposed into tangential and normal 
components; At first order, tangential deformations may always be identified 
with a reparametrization; at second order, they differ. The relationship between 
tangential deformations and reparametrizations, as well as the coupling between 
tangential and normal deformations, is examined at this order for both the metric 
and the extrinsic curvature tensors. Expressions for the expansion to second 
order in deformations of geometrical invariants constructed with these tensors are 
obtained; in particular, the expansion of the Hamiltonian to this order about an 
equilibrium is considered. Our approach applies as well to any geometrical model 
for membranes. 



PACS numbers: 87.16.Dg, 46.70.Hg 

1. Introduction 

In water, lipid molecules assemble spontaneously into vesicles which are described 
remarkably well at mesoscopic scales by a purely geometrical Hamiltonian ^|2l|3]. On 
such scales, there is a difference of several orders of magnitude between the thickness 
of the lipid bilayer and the diameter of the vesicle; it is therefore sensible to describe 
the vesicle itself as a two-dimensional surface, the relevant course grained degrees of 
freedom are purely geometrical and describe the shape of this surface. Furthermore, 
this membrane acts like a two dimensional fluid: there is no cost in energy associated 
with tangential displacements of the lipid constituents which preserve the area, and 
thus its shear modulus vanishes. In this respect, the membrane differs completely in 
its behaviour from a familiar elastic solid. As a two-dimensional fluid the membrane is 
then described by an effective energy that does not penalize tangential displacements. 
Infinitesimally, tangential displacements can be identified with a reparametrization of 
the surface. The appropriate Hamiltonian must therefore be a geometrical invariant 
under reparametrizations. 

The Helfrich-Canham Hamiltonian quadratic in the mean extrinsic curvature 
describes the penalty associated with the bending of the vesicle 0] The 
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microscopic physics of the lipid molecules is encoded in the rigidity modulus 
characterizing the stiffness of the membrane on the particular mesoscopic scale being 
considered. 

There are global constraints on the shape of the vesicle: the area is fixed, and on 
time scales relevant in experiments, the enclosed volume is also. At first order, the 
particular composition of the lipid bilayer, and in particular the asymmetry between 
the layers, is characterized by a constraint or a penalty on the total mean extrinsic 
curvature of the surface, a quantity which captures the area difference between the two 
layers. For the sake of simplicity we will restrict our attention to a minimal geometric 
model for fluid vesicles which is known as the strict bilayer couple model 0. A more 
realistic geometric model, the area-difference model, takes the bilayer composition 
into account more precisely, in particular, the difference in stretching of the individual 
layers jHlElEIj- The two models are related by a Legendre transformation; the formal 
questions we address apply to both. 

Our aim in this paper is to provide an approach towards a covariant theory of 
deformations of a membrane described by the Helfrich-Canham Hamiltonian, although 
our considerations will not depend on the details of this model; they apply equally 
well to any reparametrization invariant geometrical theory of membranes. The basic 
variables are the shape functions describing the surface. We examine how the geometry 
of this surface changes under a deformation. The approach we will adopt complements 
the one presented in where the deformation was decomposed into its tangential 
and normal components with a focus on the latter; it is simliar, at least in spirit, to 
the approach taken by Cai and Lubensky in their description of membrane dinamycs 

P21 Q3]. 

At first order in the deformation, the change in the Hamiltonian vanishes when the 
vesicle is in equilibrium. At this order, the tangential deformation of the Hamiltonian 
appears only in a boundary term, so that for the purposes of examining the equilibria of 
closed vesicles, it can always be neglected. This is consistent with our understanding 
that an infinitesimal tangential deformation is a reparametrization of the surface. 
Contrary to what one might expect, however, this identification breaks down at higher 
orders. Finite tangential deformations are not simple exponentials of infinitesimal 
reparametrizations. Nonetheless, as was shown in if one is interested only in 
fluctuations about equilibrium (so that the Euler-Lagrange equation is satisfied), 
tangential deformations remain irrelevant at second order. This is because, at this 
order, the tangential contribution to the deformation of any given term appearing 
in the Hamiltonian is proportional to the Euler-Lagrange derivative of that term. 
In equilibrium the sum of these terms vanishes: the second variation of the total 
Hamiltonian about equilibrium is thus always a quadratic in the normal deformation. 
This is the principal justification for the cavalier approach adopted in [1 1 j where 
tangential deformations are discarded from the outset; considering the effort one must 
expend to keep track of tangential contributions, the fact that no error is incurred 
represents a stroke of good luck. What it fails to do, however, even at second order is 
provide a correct expansion of individual geometrical tensors, such as the metric and 
the extrinsic curvature: tangential deformations not only contribute but also couple 
non-trivially to normal deformations and there is no justification to drop them. 

In this paper, we will examine the coupling between tangential and normal 
components at second and higher orders. We will also attempt to quantify the extent 
to which tangential deformations differ from reparametrizations. These issues do not 
appear to have been addressed before. Besides their value in point of principle, 
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there is also a practical value to understanding them: there are occasions when it 
is necessary to look beyond second order. To identify the stable deformations of a 
spherical equilibrium shape one needs to expand the Hamiltonian out to fourth order 
[141 115j in order to resolve a degeneracy occurring at second order. Helfrich and Ou- 
Yang did not attempt the full calculation focusing instead on a single mode. Whereas 
at second order they can be ignored, tangential deformations will contribute at higher 
orders to perturbations about an equilibrium, and thus they must be confronted. 
It is perhaps not too surprising that, to date, a renormalization group analysis of 
the fluid membrane model has not been attempted at two-loops. The number of 
terms involved, even in a straighforward Monge representation of deformations, is 
sufficient to discourage the feint hearted. To attempt such an exercise, with some 
hope of successfully completing it, it is important to identify underlying patterns in 
the expansion of the Hamiltonian. 

This article is organized as follows. Sect. 2 introduces the geometry of two- 
dimensional surfaces and the Helfrich-Canham Hamiltonian for lipid membranes. 
Sects. 3 and 4 consider how the intrinsic and extrinsic geometries of the surface 
change under an infinitesimal covariant deformation up to second order, respectively. 
In Sect. 5, the deformation is decomposed in its normal and tangential components and 
we examine the relationship between tangential deformations and reparametrizations 
both at first and second order. In Sect. 6 we derive the first order variation of 
the Helfrich-Canham Hamiltonian and we identify its Euler-Lagrange derivatives, 
obtaining the shape equation which determines the equilibrium configurations. In 
Sect. 7, we derive general expressions for the stresses and torques associated with the 
Hamiltonian, providing an alternative derivation of the results of [HjJ. The second 
order variation of the Hamiltonian is derived in Sect. 8. We conclude in Sect. 9 with 
some final remarks. 

2. Geometric model 

We model a lipid vesicle as a 2-dimensional surface E embedded in 3-dimensional 
space. This surface is specified locally in parametric form by three shape functions 
X = (X 1 , X 2 , X 3 ), x = X(£ a ), where the coordinates x = — (x l ,x i ,x i ) describe 
a point in space, £° = (?\£ 2 ) are arbitrary coordinates on the surface. 

First, we recall briefly some basic facts about the geometry of surfaces. For a 
thorough introduction to this subject see e.g. ffHj- The two tangent vectors 

to E are e a = 9 a X, with d a = d/d^ a . The metric induced on E by the embedding 
is defined by g a b = e a • ef,. Latin indices are lowered and raised with g a b and its 
inverse g ab , respectively. The induced metric defines the infinitesimal area element 
with dA = yfg d 2 £, where g denotes the determinant of g a b- The unit normal to the 
surface E, n, is defined implicitly by n • e Q = 0, and n • n = 1. We note that the basis 
vectors {e a ,n} are complete: given any two vectors U and V, 



The classical Gauss- Weingarten equations describe the the expansion of the surface 
gradients of the basis {e a , n} adapted to the surface E in terms of the basis: 



Here T c ab denotes the Christoffel symbols of the E covariant derivative compatible with 
g a b, such that for an arbitrary surface vector V b we have V a V b = d a V b + T b c V c . By 



U • V = (U • n)(V • n) + (U • e„)(V • e a ) . 



(1) 



d a e b = T c ab e c - K ab n , 
<9 a n = K ab g bc e c . 



(2) 
(3) 



Second variation of the H elf rich- Canham Hamiltonian. 



4 



compatible we mean V a gb c = 0. The ChristofFel symbols r^ b are given in terms of the 
induced metric by 

r a hc = e a ■ d b e c = \g ad {d b g cd + d c g bd - d d g bc ) . (4) 

Geometrically, the T ab are purely intrinsic: they depend only on the induced metric 
g ab . In a geometrical covariant description of the surface, the ChristofFel symbols 
appear only through the covariant derivative. The intrinsic Riemann curvature tensor 
of E quantifies the degree of failure of the covariant derivative V a to commute, 

(V a V b - V b V a )^ c = K c dab V d . (5) 

In terms of the Christoffel symbols, the Riemann tensor takes the form 

Tl a bcd = dcT db - d d T^ b + T^ e T e db - T^ e T e cb . (6) 

Contraction of the Riemann tensor gives the Ricci tensor lZ ab = TZ c acb , and the scalar 
curvature is given by contraction with the contravariant metric 1Z = g ab 1Z ab . For a 
two-dimensional surface, the Riemann tensor is completely determined by the scalar 
curvature 

Tlabcd = (H/2) (g ac gbd - gadgbc) , (7) 

which implies, in particular, TZ a b = (l/2)lZg a b- The scalar curvature of a two- 
dimensional surface is twice the Gaussian curvature G, i.e. 7Z= 2G. 
The extrinsic curvature tensor of the surface E is 

K ab = -n • d a e b = K ba ■ (8) 

As a real symmetric two by two matrix it can always be diagonalized. In particular, 
the eigenvalues ci, C2 of the matrix K a b = g bc K ca are the principal curvatures of the 
surface. The trace with the contravariant metric K = g ab K ab is the mean curvature 
of the surface. With respect to the principal curvatures K = c\ + c%. In the literature, 
often the mean curvature is denoted by H = (1/2)K. The Gaussian curvature is given 
in terms of the principal curvatures by their product, G = c\ C2- 

The intrinsic and extrinsic geometries of E are related by the Gauss-Codazzi- 
Mainardi equations, which arise as integrability conditions for the Gauss- Weingarten 
equations ©, 

Tlabcd K ac K bd + K ad K bc = , (9) 

V a K bc - V b K ac = . (10) 
We will use extensively their contractions with the contravariant metric g ab : 

K ab -KK ab + K ac Kb c = 0, (11) 

K-K 2 + K ab K ab = , (12) 

V b K a b -V a K =0. (13) 

Note that for a 2-dimensional surface the contracted Gauss-Codazzi equation (|12|l 
contains the same information as ©. 

The fluid state of the lipid vesicle implies that in an effective mesoscopic 
description shear is negligible, therefore the vesicle Hamiltonian has to be invariant 
under reparametrizations. Moreover, as recognized long ago the important mode of 
deformation is out of the surface, corresponding to a bending jUlSlIB]. The bending 
energy is quadratic in the mean extrinsic curvature 

F b =a I dAK 2 , (14) 
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where the constant a denotes the bending rigidity. At the same order, we have also 
the Gaussian bending energy 

Fg — a G JdAK, (15) 

with ac the Gaussian bending rigidity. However, if the surface has no boundary, by 
the Gauss-Bonnet theorem, the Gaussian bending energy is a topological invariant 
(see e.g. [TH|): 

F G = 87ra G (l-g), (16) 

where g is the genus of the surface. As such it does not contribute to the determination 
of the equilibrium configurations of the membrane. Note that at the same order we 
have also the geometrical invariant J dA K ab K a t,. However, it is not independent 
since it is related to the bending and Gaussian bending energies via the Gauss-Codazzi 
equation (jT2|) . 

The lipid vesicle is subject to various geometric constraints. The low solubility 
of the lipid molecules implies that its area A is constant. The low permeability of 
the membrane implies that the enclosed volume V is constant. We write the enclosed 
volume V as a surface integral with 

V= \ [dAn-X. (17) 



3 . 

The bilayer architecture of the lipid membrane is captured, in a first approximation, 
by a constraint on the area difference between the layers. This is expressed as constant 
total mean curvature (see jZj) 



M = J dAK, (18) 

since, as we shall see below, the difference in area is proportional to the mean extrinsic 
curvature. 

Therefore we are led to consider the Helfrich-Canham Hamiltonian 

F = F b + nA + f3M-PV, (19) 

where /j, P, (3 are the Lagrange multipliers that enforce the constraints of constant area, 
volume, and total mean extrinsic curvature, or constant area difference, respectively. 

It is important to emphasize that a more realistic model of the lipid vesicle which 
takes the bilayer architecture more accurately into account is the area-difference model 
|B1 El E3] ■ A thorough discussion of the extant curvature models for lipid membranes 
can be found e.g. in the reviews pfl lTTTl I2U] . 



3. Deformations of the intrinsic geometry 

Let us consider now deformations of the surface E. A one-parameter infinitesimal 
deformation of the shape functions X(£ a ) can be described by 

X(C) = X(C) + eW(C) ; (20) 
W(£ a ) is an arbitrary vector field, and the constant e an infinitesimal parameter. Such 
a deformation gives us a new surface S. We begin by examining how the intrinsic 
geometry of the two surfaces E and S are connected up to second order in e. A tilde 
will be used to denote the geometrical quantities that characterize E. The content of 
this section can be found in many monographs on differential geometry. However, it 
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is often presented in an abstract notation quite unfamiliar to the working physicist. 
For this reason, we offer a self-contained derivation of the relationship between the 
geometries of E and E. 

We will not consider deformations in which W itself depends explicitly on X, such 
as rigid rotations of the surface or contexts where it is useful to tie the surface local 
coordinates £ a to the embedding itself. An example of the latter is to parametrize 
the surface by arclength along some priviledged direction. Consistency would require 
then that the coordinates themselves suffer a deformation. Failure to account for this 
fact is a source of frequent errors in the literature. 

At fixed values of the arbitrary surface coordinates £ a , the tangent vectors to E 
and E are related by 

e a = e a + e a(1 ) = e + eW 0j (21) 

where we define W a = <9 a W. The change in the tangent vectors is only linear in e, 
like the shape functions X themselves. For a constant deformation, W = a = const., 
the tangent vectors coincide. A translation of the surface will not change its geometry. 
It follows, using (|21|l . that the induced metric on E takes the form 

g ab = e a ■ e b = g ab + 2e(e (a ■ W 6) ) + e 2 (W a ■ W 6 ) , (22) 

where the round brackets enclosing indices denote symmetrization, i.e. A^ ab ^ = 
(l/2)(A ao + Aba). It should be emphasized that this expression is valid to all orders 
in e; it terminates at second order. 

The area measure on E is dA = \fgd 2 ^, with g the determinant of g ab . It is 
related to the area measure on E by the expression 

Vff = + e(e a ■ W a ) + € - [(n • W a )(n • W„) 

+ (e Q • W Q )(e b • W b ) - (efc • W a )(e a • W b )]} + 0(e 3 ) . (23) 

Note that the two terms on the second line have the structure of a determinant for the 
2-dimensional matrix (e Q • W b ). In order to derive this expression, we need the inverse 
induced metric g ab , defined by g a cg bc = 5 a , together with (|2*2*|l . For this purpose, we 
expand g ab and g a b in powers of e. Collecting terms linear in e, we have the condition 
g ab (i)9bc + g ab 56c (l) = 0, where the number in parenthesis refers to the order in e. 
This gives 

g ab {l) = -2e(e< a • W 6 )). (24) 

At second order in e, we have the condition g ab (2)9bc + g ab (i)<?6c(i) + g ab gbc(2) = 0, 
which, in turn, yields 

g ab (2) = e 2 [(e a • W c )(e b • W c ) + 2(e< a • W c )(W b > • e c ) - (n • W a )(n • W b )] , (25) 

where we have used the completeness relation to get (e c ■ W a )(e c ■ W b ) = 
(W a • W h ) — (n ■ W a )(n • W b ). It follows that for the inverse induced metric on 
E we have 

~ab = g ab _ 2g ( e (a . W &) ) + £ 2[( e a . Wc )( e & . W c ) - (n . W a )(n • W & ) 

+ 2(e (a • W c )(W b) •e c )]+C(e 3 ). (26) 
Note that g ab , unlike g ao , has corrections at all orders in e. 
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We calculate 



= y/g + e 



e 2 



9V? 

<9e 2 



0(e d 



(27) 



e=0 



'9 9 



de 2 V de 
Taking the e — ■> limit, we have 

v^ (1) = e (e Q ■ W° 
At second order, we have 



[5 Qb (e a ■ W b ) + eg ab (W a ■ W b )] . (28) 



(29) 



de 2 



[T\e a ■ W b ) + e{W a • W h )] + Vff 3 afc (W a • W b 



de 



[(e a • W b ) + e(W a • W b )] . 



(30) 



Using l|2t)|l and (|28|l in this expression, and taking the e — > limit gives the second 
order correction in i|23|l . 

Finally we consider the intrinsic scalar curvature R. We restrict our attention to 
the first order correction. We use the Palatini identity for the first order correction of 
the Christoffel symbols T^ h , 

' _ ~~ ~~ (31) 



r afc(l) = 2~9 C (^b9ad(l) +V a #6d(l) -Vd9ab(l)) 



to obtain, using 122L 

T c ab{1) = e(e c - V a W b + W c - V a e b ) 

- e[e c ■ V a W b - ^ ab (W c • n)] , (32) 

where we have used the Gauss- Weingarten equation (j2J in the second line, and the 
fact that the covariant derivative is torsionless, i.e. V a W b = V b W a . We recall that 
^afc(i)' urm ke r^ b , transforms as a tensor under surface reparametrizations. Using the 
definition of the Riemann tensor given by 10, we have that at first order the Riemann 
tensor is 

Tl a bcd(i) = V c r^ b (i) - V d r° b(1 ) , (33) 
so that, inserting (|32|l . we have 

n a bod(l) =€[- Tl e bcd {e a ■ W e ) + TZ cdb e (e e ■ W a ) + K bc (n ■ W d W a ) 

- K c a (n ■ V d W b ) + K d a (n ■ V c W b ) - K db (n ■ V C W Q )] , (34) 

where we have used the definition of the Riemann tensor and both the Codazzi- 
Mainardi equations (|10|l and the the Gauss-Codazzi equations ®. For the Ricci 
tensor this implies 

Kbd(i) = V c rg d ( 1 ) - v d r^ b(1 ) 

= 4-(K e bad + n e dab )(e a ■ W e ) + K ab (n ■ X7 d W a ) 

+ K a a (n ■ V a W b ) - K db (n ■ V C W C ) - K(n ■ V b W d ) , (35) 
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These expressions are valid for a hypersurface of arbitrary dimension. Restricting our 
attention to 2-dimensional surfaces, and exploiting the fact that both the Riemann 
tensor and the Ricci tensor can be expressed in terms of the scalar curvature 1Z, we 
have that 

9 ab n ab{1) = e[2(K ab - Kg ab )(n ■ V Q W b ) - TZ(e a ■ W a )} , (36) 
For the scalar curvature we have then 

= g ab (i)T^ab + g at "R-ab(i) 

= 2e[(K ab - Kg ab )(n • V a W fc ) - K(e a ■ W«)] . (37) 

Note that it depends on two derivatives of the deformation vector W, and that it 
involves the extrinsic geometry of the surface. 

In general, it is always possible to expand any geometrical quantity / on the 
deformed surface £ in terms of W. With / = /+/(i)+/(2)j besides direct computation, 
we can obtain the second order term /( 2 ) by 'deforming' the first order term /m, that 
is via the important identity 

/(2) = ^ /(l) ■ ( 38 ) 

Here /m is to be understood as the expansion to order e of /m, as we illustrate below. 
This alternative approach is particularly useful when we consider global geometric 
quantities associated with E (see Sections 6,8). Note that, in agreement with the 
identity 1)38(1 . we have, for example, 

y/9(2) = \7~9ix) = I [V? (1) (e a • W ) + y/gg$(e a ■ W b ) + Vff(W a • W a ) 



y (i) 

V5[(e Q • W a ) 2 - (e a • W b )(e a ■ W 6 ) - (e 6 ■ W a )(e a ■ W b ) 



2 

+ (W Q • W a )] , (39) 

which, using the completeness relation in the last term, reproduces the second order 
contribution to l|23|l . The identity l|38|l can be proved using variational techniques. 
It does not appear to be available in the literature. Of course, this could be just a 
shortcoming of our search. 



4. Deformation of the extrinsic geometry 

Let us turn now to the extrinsic geometry of the deformed surface S. For its unit 
normal n, we use the defining relations n • e a = 0, n • n = 1, together with (|21|l . We 
expand n and we obtain the relations 

e(n- W a ) + (n (1) ■ e a ) =0, 

e(n ( i) • W a ) + (n (2) ■ e a ) = 0, 

n-n (1) =0, 

n(i) • n(i) + 2n • n ( 2) =0, 

which provide six equations for the six unknowns n^) , n( 2 ) . Some simple algebra gives 

n = n e(n ■ W a )e a + e 2 [(n ■ W fc )(e b • W Q )e a - ±(n ■ W a )(n ■ W„)n] + 0(e 3 ) . (40) 

Note that if the deformation is such that (n ■ W a ) = the normals to the two surfaces 
coincide. This happens for parallel surfaces, defined by W = an, with a constant pi]. 
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The extrinsic curvature of S is K ab := — n • <9 a e b . Expanding the right hand side 
to second order in e, we obtain 

K ab = K ab - n (1) • d a e b - n • d a e b (1) - n (1) • <9 a e (1) - n (2) ■ d a e b (1) + 0(e 3 ) . (41) 
We use ipil) . (14011 . together with the Gauss- Weingarten equation for E, J5J). We obtain 

iL 6 = ^afc-e(n-V a W b ) + e 2 [(n-W c )(e c -V a W^ 

We note that K ab transforms covariantly under reparametrizations of £ and that it 
involves two derivatives of the deformation vector W. Note that for parallel surfaces, 
one has that K ab ^ = eaK a c K cb and K ab ^ 2 ) = 0. 

For the trace of the extrinsic curvature K, Eqs $261 . (|42[1 imply 

K = K - e[(n • V a W a ) + 2K ab (e a ■ W b )} 

+ e 2 [2(e a ■ W b )(n • V Q W b ) + 2K ab (e a ■ W c )(e c • W b ) 
+ K ab (e a ■ W c )(e b • W c ) + (n • W c )(e c • V a W a ) 

- K ab (n ■ W a )(n . W 6 ) - y (n ■ W c )(n • W c )] + 0(e 3 ) . (43) 

At this point we can use the Gauss-Codazzi equation l|12fl to check the validity 
of these expressions for the deformation of the extrinsic curvature. At first order, 
using i|26|l . I|42|l . (|43|l . we reproduce l|37|) . Note that the check is non trivial; 
it requires extensive use, in various degrees of contraction, of the Gauss-Codazzi- 
Mainardi equations. 



5. Deformations of the geometry: decomposed 

The existence of a basis adapted to the surface S, {e a ,n}, suggests a natural 
decomposition of the deformation vector W into its tangential and normal components 
(see e.g. [Ill), 

W = $ Q e a + $n, (44) 

It follows that the projections of the first two derivatives of the deformation vector 
are 

e„ • W b = V b $ a + <PK ab , (45) 
n W a = V a <S>-$ b K ab) (46) 
e c • V a W b = V a V b $ c - <S> d K bd K ac + 2iv c(Q V h) $ + $V a if bc , (47) 
n • V a W b = V Q V b $ - <S>K ac K c b - 2K c{a V b) <S> c - <$> c V c K ab , (48) 

where we have used the Gauss- Weingarten equations J3J), ©, an d the contracted 
Codazzi-Mainardi equation (|13|l . 

Using these expressions, the basic geometric quantities that characterize £ take 
the form, to first order in e, 



9ab(l) 


= e (2K ab $ + 2V (a $ b) ) , 


(49) 




= eV5 (i^$ + V Q $ a ), 


(50) 


Kab(l) 


- e (-V a V b $ + K ac K c b t> + $ c V c X ab + 2^ c(a V b) $ c ) . 


, (51) 




= e (-V 2 $ - K ab K ab >$> + <P C V C K) , 


(52) 




= e [2(K ab - Kg ab )V a V b <5> + $ c V c ft] , 


(53) 
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We can thus identify both a normal (linear in $), and a tangential deformation (linear 
in <& a ) of these geometrical quantities. We note that these expressions coincide with 
the ones obtained e.g. in 11 . 

In particular, to first order the tangential components correspond to an 
infinitesimal (active) reparametrization of the surface. Indeed, each of the three surface 
scalars X(£ a ) transforms as 

6 rep X = v a (Z b )d a X = v a e a , (54) 

under a reparametrization £ a — ► £ a — v a . This is exactly the effect of a tangential 
deformation at first order with the identification of v a with the surface vector field 
defined by the projection $ a = W ■ e°. Reflecting this fact, at first order, geometrical 
quantities transform as surface Lie derivative along the surface vector field <I> a . For 
example, setting $ = in i|49|) . we have that 

5ah(i)tang. = 2V (a $ h) = C^g ab , (55) 

where £$a denotes the surface Lie derivative along $ a . 

It is important to emphasize that, having made this identification at first order, 
the tangential part of any total geometrical invariant of £ is always given by the 
integral of a total divergence, which vanishes over a closed surface without boundaries. 
To see this, consider an invariant / = J eL4./(£°), where f(£, a ) is a scalar under 
reparametrizations. Since / is a scalar we have simply that 

/(i)tang. = ® a d a f, (56) 
no matter how complicated the dependence of / on the geometry might be. Moreover, 
setting $ = in (5Dt , we have v^-^tang = v 7 ^ ^a$ a , and therefore 

*(i)tang. = / dAV a (<S> a f). (57) 

At first order, we can always disentangle the physical normal deformation from 
reparametrizations and we can safely set $ a to vanish. Matters, however, are not 
so simple at second order. Let us consider the second order variation of the metric. 
Using the completeness relation Q it takes the form 

9ab (2) = e 2 [(e • W b ) (e a • W b ) + (n • W a ) (n • W a )} 

= e 2 [(Vfc$ Q + $K afc )(V b $ a + <S>K ab ) + (V a $ - $ c ^ ac )(V a $ - $ d K ad )] 

= e 2 [V Q $V b $ + K a c K cb $ 2 + 2^ c(Q V b) ($$ c ) 

+ V a $ c V 6 $ c + K ac K bd $ c $ d ] , (58) 

and with (|49|l this completely describes the deformation of the metric to all orders. 
At second order, normal and tangential deformations begin to talk to each other. 
When both normal and tangential deformations are present, there is a mixing. The 
purely tangential deformation at this order is certainly not simply a second order 
reparametrization, in the sense of a composition of Lie derivatives. This might appear 
to be obvious: the second order tangential deformation involves the extrinsic geometry 
through the quadratic term in K ab , whereas a reparametrization is a purely intrinsic 
concept and as such it should not involve K ab . One has to be careful, however, to 
check how the dependence on the extrinsic geometry enters. 

Let us look more closely at the issue of reparametrization covariance at second 
order. Consider an infinitesimal change of coordinates on the surface £ a — > £' a — 
£ a + v a , and a surface scalar field f(£, a ). By definition of scalar field we have 



Second variation of the H elf rich- Canham Hamiltonian. 



11 



/'(£' a ) = /(£")■ I n order to evaluate the change in the scalar field at the same 
point, we expand /'(£'") = /'(£ a + v a ) in powers of v a , 

f(e a ) = rm + v a d a f(o + l(v a d a )(v b d b )f(n +■■■, m 

where we replace /' by / in the last term since it is already of second order in v a . For 
the middle term note that 

v a d a f(c) = v a d a {f(e a ~v a )} 

= v a d a [f(e)-v b d b f{e)] 

= v a d a [f{C)~v b d b f{C)\- (60) 
Therefore, to second order, we have 

&ep/(0 = f'(C) - f(C) = -v a d a f(C) + \{v a d a ){v b d b )f{e) ■ (61) 

At first order, we have minus the Lie derivative of the scalar field, since now we are 
considering passive transformations. The second order contribution is the composition 
of two Lie derivatives. 

In particular, for the embedding functions we obtain 

JrepX = -v a e a + 1 (v a V a v»e b - v a v b K ab n) . (62) 

The important point here is that, at second order, a reparametrization will generally 
alter the embedding functions. In contrast, by construction X is only modified at first 
order in W. Note also that, at second order, a reparametrization produces a change 
in X along the normal direction. Moreover, it depends explicitly on the extrinsic 
geometry. This justifies our earlier caveat. 

The tangent vectors e a transform as covariant surface vectors under 
reparametrization: we have at first order for a vector f a 

SrepL (i) = -v b d b f a - d a v b f b . (63) 

Thus for e a : 

^repe a (l) = - v b d b e a - d a v b e b 
= - v b V b e a - W a v b e b 

= v b K ab n - V a v b e b , (64) 

involving both tangential and normal parts. Note how the normal contribution 
projects out of <5 re pffafc(i) = 2e (a ■ <5 r epe 6 ) (i) = 2V (a u 6) . At second order, 

^repe Q ( 2 ) = - (-v b d b S T epe a (1) - 5 T epe b (i)d a v b ) . (65) 
Similarly, we obtain for the metric at second order 

£repff a &(2) = ^K(V c V a v fc + V c V b w a ) + (V c u a )(V b w c ) 

+ (VcV b )(V a v c ) + 2{V a v c )(V b v c ) . (66) 

Note that <5rep5'ab(2) is manifestly intrinsic and, as such, distinct from the contribution 
t° 5ab(2) quadratic in $ a on setting $ a = v a . We note also that this expression 
coincides with the metric at second order induced by (|62|l 

^rep5a&(2) = 2e (a • ^repef,) (2 ) + <5 r epe a(1) • c>rep e b(i) ■ (67) 
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At second order, we see that we cannot disentangle the physical normal deformation 
from reparametrizations and we cannot set $ a to vanish. However, as we will see 
below in Sect. 8, when considering the second order deformation of global geometric 
invariants, the tangential component of the deformation will appear only in boundary 
terms or in terms that vanish when the membrane is at equilibrium. 



6. Variation of the Helfrich- Canham Hamiltonian: first order 

Let us expand the Hamiltonian F[X] as given by H19|) to first order in e. This will 
allow us to identify its Euler-Lagrange derivative, and the equilibrium conditions 
for the vesicle. (For a different approach emphasizing the normal component of the 
deformation, see 

The expansion can always be written in the form 

F (1) [X,W] = e J dAE F ■ W + e J dAV a Q a . (68) 

Here Ej? denotes the Euler-Lagrange derivative for -F[X]. The quantity Q a appearing 
in the total divergence in the second term is the Noether current |16| . which will be 
used below in Sect. 7 to derive the stresses and torques acting on the surface associated 
with F[X\. 

We use the results of Sect. 2 to derive, term by term, the various contributions 
to H68|) . For the area of the vesicle, we find using l|23|l that 

A (1) =eJdA(e a -W a ). (69) 

To cast this expression in the form l|68(l . we integrate by parts, and obtain 

A (1) = e J dAKn-W + e J dA V a (e a -W). (70) 

Therefore, the Euler-Lagrange derivative of the area is purely normal, and proportional 
to the mean extrinsic curvature, E^ = = Kn. The first feature is common to 
all reparametrization invariants: to first order in e, tangential deformations contribute 
only boundary terms, as shown in the previous section. The latter tells us that minimal 
surfaces, extrcmizing the area, have vanishing mean extrinsic curvature, Ea — K — 0. 
Note that for a constant normal displacement W = an we have 

A (1) =a J dAK = aM , (71) 

so that the total mean extrinsic curvature is proportional to the area difference in the 
normal direction. 

If we require that the area be infinitesimally locally invariant, then we have a 
constraint on W of the form e a • W a = 0. This does not, however, alter the value of 
the Euler-Lagrange derivative. 

For the volume enclosed by the vesicle, we use the definition Q17[l. together with 
ijltl jl . to derive 



e 

= 3 



J dA {1) (n ■ X) + J dA (n (1) • X + en • W) 
dA [(W • n) + (e a • W Q )(n • X) - (n • W a )(e a • X) 
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We integrate by parts the second and third terms and neglect a total divergence to 
obtain 

V ( i) = e J dA n ■ W , (72) 

therefore we find that the Euler-Lagrange derivative of the enclosed volume functional 
is simply unity, Ey = 1. 

Let us consider now the total mean extrinsic curvature, M, as defined in 118|) . 
We use (GUI and g3J to derive 

M (1) = e J dA [K(e a ■ W«) - 2K ab (e a ■ W b ) - (n • V Q W a )] . 

To put it in the form ((68(1 . we integrate by parts and use the Gauss- Weingarten 
equations J2J), © to arrive at 



M (1 



(i) 



e J dA Tl (n ■ W) + ej dA X7 a {(Kg ab - K ab )(e b ■ W) - (n • W a )] . (73) 

The scalar intrinsic curvature appears as the Euler-Lagrange derivative of the total 
mean extrinsic curvature functional Em = TZ-- 

As mentioned above, the Gaussian bending energy F G as given by l|15fl is a 
topological invariant. As such, we expect both Fg(\) an( l ^G(2) to be a total 
divergence. As this provides a non-trivial check, let us consider its expansion to 
first order. Moreover, in any case, we are interested in the non-vanishing contribution 
to the Noether charge. Using Ij98(l . Ij37(l . we have 

F G (1) = ea G J dA [2{K ab - Kg ab )(n ■ V a W b ) - TZ(e a ■ W„)] . (74) 

We integrate the first two terms by parts, and use the the second Gauss- Weingarten 
equation J3Jl, together with the contracted Codazzi-Mainardi equations 1)13(1 . to obtain 

F G (1) = ea G J dA [-TZg ab + 2KK ab - 2K ac K b c ](e a ■ W b ) 

+ ea G J dAV a [2(K ab - Kg ab )(n ■ W b ) - TZ{e a ■ W)] , (75) 

where the first line vanishes because of the contracted Gauss-Codazzi equation i|ll|) . 
Therefore F G ( X ) is given by a total divergence. With the help of the relationship 1(38(1 . 
the second order term F G ( 2 ) is a total divergence as well. To see this, recall the fact 
that the 'deformation' of the divergence of a vector density is equal to the divergence 
of the deformation. In this particular case we have 

F G{1) = ea G J dA V a Q a G = ea G J d 2 £ d a (^ Q G ) , (76) 

where Q G denotes the contribution of the Gaussian bending rigidity to the Noether 
current, and it is given explicitly by the argument of the covariant derivative in 1(75(1 . 
At second order we have then the total divergence 

F G (2) = \f G {x) - \c*g J d 2 £ d a [{Jg Q G ) {1) ] . (77) 

Finally, for the bending energy 1(14(1 . we obtain, using 1231) . 1(431) . 

F H1) = efdA [K 2 (e a • W a ) - 4KK ab (e a -W b )-2K(n- W a W a )} , (78) 
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and integration by parts twice gives 

F b (1) = e J dA [-2V 2 K + K 3 - 2KK ab K ab ] n • W 

+ e J dAV a [(K 2 g ab - 2KK ab )(e b • W) + 2(V a K)(n • W) - 2K(n ■ W a )] (79) 

Therefore the Euler-Lagrange derivative of the bending energy functional is 

E Fb = - 2V 2 K + K 3 - 2KK ab K ab 

= -2V 2 K + K(2K-K 2 ), (80) 

where we have used the Gauss-Codazzi equation ()12(l to obtain the second line. 

We are now in a position to write down the equilibrium conditions for the 
Hamiltonian JT^J. We set 

E F = E F n, (81) 

where 

E F = a [-2V 2 K + K(2K - K 2 )] + fj,K + (3K - P . (82) 

Then equilibrium configurations that extremize the Helfrich-Canham Hamiltonian 
satisfy 

E F = 0. (83) 

This is known as the shape equation . Note that it is a nonlinear 4th order partial 
differential equation. Progress in the understanding of its space of solutions has been 
limited to the special case of axisymmetric configurations, see e.g. |2U| . 

7. Stresses and torques 

The Noether current appearing in 1)68(1 allows the derivation of the stresses and the 
torques acting on the membrane, using the invariance under rigid motions in space. 
This was done in jT^J by decomposing the deformation in its normal and tangential 
parts. Here, we provide an alternative derivation, which has the advantage of being 
more direct. See also |2"21l%3"] . 

Before, we proceed, for the purposes of this section it is convenient to rewrite the 
Helfrich-Canham Hamiltonian (|19|l by isolating the volume term as 

F = F S -PV , (84) 

with the surface part of the Hamiltonian F s = F b + Fq + [iA + (3M . 

Collecting the various total surface divergences appearing in (|T0|) , (f?3"|l , lfT5|l , lt7^|) , 

it is straightforward to identify the Noether current associated with F s , 

Q a = a [{K 2 g ab - 2KK ab )(e b ■ W) + 2{V a K){n ■ W) - 2K{n ■ W Q )] 
+ f3[(Kg ab - K ab )(e b • W) - n • W a ] + /i (e a • W) 

+ 2a G (K ab -Kg ab )(n- W 6 ) . (85) 

Note that although the Gaussian bending energy does not contribute to the Euler- 
Lagrange equations, it does appear in the Noether current <2°. 

We consider now a (simply connected) piece of the membrane, which we denote 
by So, bounded by a curve C, and we specialize the variation of the Hamiltonian (|68|l 
to this arbitrary region of the membrane. We have that 

F s (1) = e f dA [ E F , n • W + V ' a Q a ] . (86) 
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We exploit the invariance of the Hamiltonian under rigid motions in space. First, we 
consider an infinitesimal translation W = a, with a constant. As the Hamiltonian is 
invariant under translations, the left hand side of l|86l) vanishes, and with the stress 
tensor f a defined by 

Q a = a • f a , (87) 
it follows that we can write the Euler-Lagrange derivative as a conservation law 

E Fb n = V a f a , (88) 

where the stresses associated with the Hamiltonian (|19|) are given by 

f a = -a [{K 2 g ab ~ 2KK ab )e b + 2(V a X)n] - (3{Kg ab - K ab )e b - n e a . (89) 

We emphasize that it is far from obvious from the shape equation itself (|83|l that it 
can be written as a conservation law. 

There are three conservation laws, and only one shape equation. This is a 
consequence of the reparametrization invariance of the Hamiltonian. This statement 
can be made explicit using the decomposition of the stress tensor into tangential and 
normal parts as follows: 

f a = / Qb e b + rn. (90) 

The surface covariant derivative then gives 

V.f -K ab f ab =E Fs =P, (91) 

V a f ab + K b a f a = 0. (92) 

The first equation is the shape equation expressed in terms of the projections f a and 
f ab . The second equation expresses the content of reparametrization invariance as a 
consistency check: the normal stress and the tangential stress must balance exactly in 
this way. Note that this identity is potentially useful in numerical simulations, where 
reparametrization invariance is necessarily lost, and one is interested in quantifying 
the degree of violation. 

The physical meaning of the stress tensor f a is perhaps best illustrated by 
considering the total force per unit length acting on the curve C. Concretely, C 
may be the shape of an edge of the membrane |24| . or the line boundary between 
the two phases of a two-component vesicle |2S1- If we consider a basis {t,l} on the 
surface adapted to the curve C that bounds So 5 with t tangent to C, and 1 = l a e a the 
(outward) normal to C on the surface, we obtain the force per unit length acting on 
C, /(jf = f , as 

f = [K [l± (2aK + 0)] t + [aK(K±_ - K\\) - f3K\\ - fx] 1 - 2a(V_LiT)n , (93) 

where we denote the projections of the extrinsic curvature onto the surface as 
K\\ = K ab t a t b , K x = K ab l a l b , and K ± \\ = K ab t a l b . Note that K = K x + K\\, 
and 1Z — 2(K\\K± — K±n)- = l a V ± denotes the covariant derivative along the 
direction normal to the curve C. 

Similarly as we showed for translations, for an infinitesimal rotation of the form 
W = b x X, we can obtain the torques acting on the surface associated with the 
Hamiltonian Ijl9|l . We define the total angular momentum m a 

Q a = -b • m a , (94) 

where the torque m a can be split in its 'orbital' and 'differential' parts as 

m" = Xxr + s°. (95) 
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From the Noether charge (|85|l we obtain directly that the terms contributing to s a are 
the ones involving derivatives of the deformation vector, W a , so that 

s a = [(2aK + (3)g ab + 2a G {Kg ab - K ab )]e b x n . (96) 

Note that it is tangential to the surface. 

The differential torque and the stress tensor are related by \Fg\ 

V a s a = f a x e a . (97) 

We emphasize that this expression is valid also when not in equilibrium. 

8. Variation of the Helfrich-Canham Hamiltonian: second order 

In this section, we exploit the general results of Sects. 2, 3 to derive the expansion 
to second order in e of the Hamiltonian 1)19(1 . As we did at first order, we derive the 
various terms that contribute to it. There are two possible strategies: on one hand 
we can perform a direct expansion, alternatively we can exploit the identity 138(1 . 
and deform the first order terms we have obtained in Sect. 6. We will adopt the 
most convenient strategy for each term. This part is a straightforward calculation. 
However, both in order to have a better understanding of the final result and to make 
contact with the variational approach of |llj . we decompose the deformation vector 
W in components. This is a straightforward calculation as well, but it turns out that 
it is possible to organize the result in a way which isolates boundary terms and a 
contribution proportional to the Euler-Lagrange derivative of each term. 
For the area, we have immediately, using (J33J), that 

A {2) = € - fdA[(n- W Q )(n • W a ) + (e a • W a )(e h • W b ) - (e 6 • W a )(e a • W b )] . (98) 

The direct expansion at second order of the volume is quite complicated. It is 
preferable to expand the first order term, so that, using the identity ((38(1 and l(72|) . we 
have 

V(2) = xV(i) = ~ 



J dA C i)(W • n) + J dA(W • n (1) ) 

= L J dA[(e a • W a )(W • n) - (e„ • W)(W a • n)] . (99) 

Let us consider the second order term in the expansion of the total mean curvature 
M. We have that 

M m = j [dA {2) K + dA {x) K {x) + dA K {2) ] , (100) 
so that using ((231) and l(43|) . we obtain 

M (2) = e 2 J dA [~(e a • W a ) 2 K - ~(e a • W b )K(e b ■ W a ) - (e a ■ W a )(n • V a W a ) 

- 2^ fcc (e a • W a )(e h • W c ) + 2(e a • W fc )(n • V a W t ) + 2^ ab (e Q • W c )(e c • W 6 ) 

+ K ab (e a ■ W c )(e b • W c ) + (n • W c )(e c • V a W a ) - K ab (n ■ W a )(n • W b )] .(101) 

There is no obvious simplification. On the other hand, using the identity ((38(1 and 
1(73(1 . we obtain, up to a total divergence, the simpler expression 



M( 2 ) = f-~- 



M (1) = jJ dA i 2 ( Kah - K 9 ab )(n ■ V a W b )(n • W) 

- K[(n ■ W a )(e a • W) + (n ■ W)(e a ■ W a )]} . (102) 
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These two expressions differ only by a total divergence. However, it is quite involved 
to extract it from (|1(J1[I . 

For the bending energy, we have that a direct expansion is the more convenient 
approach and with 

F b (2) = J [dA {2) K 2 + 2dA (1) K K (1) + dAK {1) K {x) + 2dA KK (2) ] , (103) 
using iffify . IpH^ . we obtain 



F H2) =e 2 



J dA [(n • V a W a ) 2 + 2{2K ab - Kg ab ){e a • W b )(n • V C W C ) 



+ AK(e a • W b )(n • V a W b ) + 2K(n • W c )(e c • V a W a ) 

- 2KK ab (n • W„)(n • W b ) - i^ 2 (n • W«)(n • W a ) 

+ (4K ab K cd + 4KK ad g bc - 4KK cd g ab + 2KK ac g bd + ^K 2 g ab g cd 

- \K 2 g ad g bc )(e a ■ W b )(e c • W d )] . (104) 

These expressions provide directly the second variation of the Helfrich-Canham 
Hamiltonian in terms of the deformation vector W and its first and second derivatives. 
It is desirable, however, in order to make contact with the expressions derived in jllj . 
to decompose W into tangential and normal. This involves plugging l|45|) to (|48H into 
the covariant expressions we have derived. For example, for the area functional we 
obtain 

A m = - 



J (L4[V a $V a $ + (K 2 - K ab K ab )§ 2 - 2/sr afc V ( $f,$) + K ac K b c $ a <$> b 

+ 2 J ftT$V a $ a + (V a $ Q ) 2 - (V a $ b )V 6 $ a - 2!f afc $V a $ 6 ] . (105) 

This is not the most useful form, however. Integrating by parts and isolating a total 
divergence, we can write it down in an alternative way as 

A {2 ) = y J cL4{-$V 2 $ - K ab K ab <S> 2 
+ K(K<I> 2 + K ab $ a <$> b - 2$ Q V a $) 

+ V a [$V a $ + 2$$ b (Kg ab - K ab ) - $ b V 6 $ a + $ a V b $ b ]} . (106) 

The first line is the normal part of the second deformation. The second line is 
proportional to the trace of the extrinsic curvature K. It is essential to recognize that 
K is the Euler-Lagrangian derivative for the area functional, Ea = K. Therefore, at 
equilibrium, the second line will vanish. The third line is a total divergence and, for 
a closed vesicle without boundaries, it can be set to vanish. 

This example indicates that it is possible to obtain in a systematic way simpler 
expressions by isolating terms that are total divergences. Let us use the expression (1681 
for the first variation of the Helfrich-Canham Hamiltonian (see ^T] for an equivalent 
argument in an alternative language). Using the identity (|38|l . we have that 

F (2) [X, W] = | [X, W] = |{ J dA E F (1) (n • W) + J dA E F (n (1) • W) 

+ j dA (1) E F (n-W)+J d 2 £ a Q [(V5Q a ) ( i)]} , (107) 
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where we have rewritten the second term of <|68[) so that ^fgQ. a is a scalar density 
of weight one; its divergence is then independent of the afhne connection, so that 
variation and derivation commute. Ep m denotes the first order variation of the 
Euler-Lagrange derivative appearing in (|83|l . We now use 1|23|) . I|4U|I to obtain 

F (2) [X,W] = y j dAE F [(e Q -W a )(n.W)-(e Q -W)(n-W a )] 

+ iJdAEp (1) (n-W)+ e -J d^d a [(^Q a )(i)}. (108) 

We now use the results of Sect. 4, to express in components the deformations. 
Furthermore we split the first order correction of the Euler-Lagrange derivative Ep 
in its normal and tangential parts as E F ^ — Ep ^ perp. + -^(ljtang ; where, since 
Ep is a scalar, £7? (l) tang = & a V a Ep. It follows that we can rewrite (|108fl as 

F P)[ X,W] -tjtAB, \K*> + *V„*» - *°V„* + K***\ + **°V„E F 

+ ||«E, ( „pe,p(ii-W) + i| <Z ! £ S„[(V?2")<1)] , (103) 

or as 

F {2) [X, W] = ~fdAE F (1) perp (n • W) + L J dA E F [K<£ 2 - 2$ a V a $ + K ab <S> a <S> b ] 



c 2 



2 J dA V a (E F ^ a ) + | d a [(^Q a ) (1) ] . (110) 

The second line is a total divergence, and it can be set to vanish safely. In the first line, 
the term proportional to the Euler-Lagrange derivative Ep is surprisingly simple. We 
recognize the same structure that appears in the second variation of the area functional 
in the form iJTDSJl. 

When the shape equation Ep = is satisfied, up to a total divergence, we have 
that the second variation is simply 

F {2) [X, W] = | / 'dAEp (1) pe rp . (n • W) . (Ill) 

If we set the tangential part of the deformation to vanish, <& a = 0, then the the 
second variation takes the form 

F {2) [X, W]= 6 -JdA [E F (1) perp. (n • W) + eE F K(n ■ W) 2 ] . (112) 

Let us consider now the expressions in components of the remaining terms in 
the second variation of the Helfrich-Canham Hamiltonian in the form (|110fl . For the 
volume we have that 

e 2 r 

V ( 2) = - dA [K(K® 2 + K ab $ a <S> b - 2$ a V a $) + V a ($$ a )] . (113) 

This case is quite special, since the contribution of the volume to the Euler-Lagrange 
derivative is Ey = 1, therefore, for the volume, E v m = 0. 

For the total mean extrinsic curvature specializing to components either or 
Ijl02|l . we obtain, up to a total divergence, 

M (2) = e 2 J dA [{K ab - Kg ab )§V a \7 b § - \^RK<$> 2 

e 2 f 

+ — dA K(K$ 2 + K ab <f> a $ b - 2$ Q V a $) . (114) 
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Note that the second line is proportional to the Euler-Lagrange derivative for M, since 



For the bending energy, a direct specialization to components of (|104|) produces 
74 terms, and it is impossible to tell the trees from the forest. However, the 
general considerations that lead to l|110[) imply that the dependence on the tangential 
component of the deformation $ a is determined. Therefore, we can keep only the 
normal part of the deformation so that, up to a total divergence, we obtain 



This corresponds to the form H112|) of the second variation. 
9. Concluding remarks 

We have presented in some detail a fully covariant approach to the deformations of the 
Helfrich-Canham Hamiltonian; where applicable, we have compared it to the approach 
adopted in ^I] in which tangential deformations were discarded. 

On balance, we feel that there is something to be learnt from both approaches; the 
reader who has considered both may better judge which form of perturbation theory 
is more appropriate to the issue being addressed. For example, whereas it is trivially 
obvious in the covariant approach that the metric tensor is subject to variations 
of second order and no higher, once the decomposition has been effected, this fact 
becomes heavily disguised and would appear to involve miraculous cancellations. On 
the other hand, the second order variation of the Hamiltonian about an equilibrium 
configuration is considerably more transparent when expressed in terms of normal 
deformations. 

The nature of tangential deformations of individual geometrical tensors has been 
clarified at second order, an issue clearly beyond the scope of the analysis in 
to address. We have shown giving explicit examples that, at this order and higher, 
tangential deformations are not reparametrizations. We suspect that there is much 
still to be pinned down on the issue. Even in the well explored field of general relativity, 
disentangling coordinate artifacts from physical perturbations at second order remains 
a vexed issue |26) . 

We have not examined explicitly any order in perturbation theory higher than 
second. This is not going to be simple; it remains a highly non-trivial challenge. The 
systematic approach we have outlined will, it is hoped, provide a few reliable signposts. 
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